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АВ5ТКАСТ 


A mathematical model for the simulation of non-catalytic 
solid-gas reactions in a fluidized bed reactor is proposed. 
The performance of the model in predicting solid reactant 
conversions for an ore roasting process is investigated 
using available literature data. Model development required 
simplifying assumptions. The sensitivity of the model to 
certain of these assumptions is investigated. 

Comments on the adaptability of the model for use in 
the design and study of a fluidized bed shipboard waste dis- 


posal system are made. 
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I. INTRODUCTION 


In the late sixties, the Navy began a study of shipboard 
waste disposal systems which could be used to combat a grow- 
ing pollution problem. One proposal considered a process 
for the combustion of solid waste material in a fluidized 
bed reactor. This proposal was opposed by some naval 
designers on the grounds that sufficient knowledge of the 
Operational characteristics of fluid bed reactors was not 
available, thus prohibiting accurate design and control of 
such a system. | 

This lack of fundamental design knowledge is a problem 
which has plagued fluid bed technologists for some time. 
Although the fluid bed reactor finds broad usage in the 
chemical industry today, its complex nature makes accurate 
mathematical modeling difficult. Design and control of 
fluid bed systems in the past has, therefore, been based on 
the application of operational engineering experience. 

This work was undertaken with the goal of studying and 
developing a mathematical model of a fluidized bed reactor. 
The Bubble Assemblage Model proposed by Wen and Yoshida [1] 
Provided the starting point for this investigation. 

If a reasonable model could be developed, it could be 
used as a design tool and in studies of the operational 


characteristics of the reactor system. 
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II. LITERATURE SEARCH 


Early tracer studies [2] on large scale fluid bed reac- 
tors revealed that tracer concentrations in the exit gas 
stream exceeded concentrations in samples taken directly 
from the bed. This observation led to the development of 
the two phase model of fluidization which predicts that the 
major portion of the gas passes through the bed in the form 
of bubbles. The.primary developers of this theory include 
Toomey and Johnstone [3], Shen and Johnstone [4], Pansing 
[5], and Lewis et al [6]. 

The two phase model pictures the bed as consisting of 
two single phase reactors in parallel as shown in Figure 1. 
imc o1, the lean or bubble phase, represents the primary 
means of gas throughput. It is generally treated as being 
in plug flow and containing no solids. Phase 2, the dense 
or emulsion phase, contains the solids surrounded by inter- 
Stitial gas. Assumptions on the mixing patterns in the 
emulsion phase range depending on the investigator. Mathias 
and Watson [7], Massimillia and Johnstone [8], and Gomez- 
plata and Shuster [9] assumed plug flow behavior while Shen 
and Johnstone [4] and Lewis et al [6] have studied condi- 
tions of plug flow and complete mixing. The model is also 
characterized by a gas exchange coefficient which accounts 
for gas flow between the phases. 

May [10] extensively studied the classic two phase model 


and was the first to introduce diffusion type coefficients 


T5 





Emulsion 
Phase 


Bubble 
Phase 





Figure 1. Six Parameter Two Phase Model. 
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to characterize the vigorous mixing in the emulsion phase. 
This approach was investigated further by van Deemter [11] 
who applied the model in a study of gas mixing experiments 
in order to determine the diffusion and gas exchange 
@eerfticients. 

The mathematical formulation of the two phase model leads 
to two partial differential equations containing six para- 
meters and representing the system material balances. A 


general expression of the model equations can be given as 











2 
2008 97008 2068 _ 
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Their solution has been a center of controversy because of 
different approaches to the choice of boundary conditions. 
McCracken [12] obtained solutions for various mixing 
patterns and in a second work [13] reviews other numerical 
solutions of the model equations. 

The two phase model, although representing a useful 
engineering approach to the understanding of fluid bed 
reactors, suffers from serious disadvantages. The model 
parameters are, in a sense, adjustable and, therefore, while 
fitting experimental data well for particular cases, their 
generality for use in scale up and design are questionable. 
Furthermore, the description of a uniform bubble phase and 
constant gas exchange coefficient ignores effects of chang- 


ing bubble size on bed operation. 
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Pioneering work in development of a more general model 
was completed by Davidson and Harrison [14] in 1962. These 
investigators theorized that fluid beds could be accurately 
described by the application of the fundamental properties 
of bubbles rising in a fluidized medium. Their work and a 
recent text by Kunii and Levenspiel [15], presents a compre- 
hensive survey of the present state of knowledge in this 
area. 

The second generation models, in general, represent the 
fluid bed as a three phase system. The dense phase, is split 
to include a cloud phase which represents solids directly 
influenced by rising bubbles. Entrapment of solids in the 
wakes of rising bubbles is also considered. Information 
on bubble size and rise velocity, cloud development, gas 
interchange rates etc., derived from independent investiga- 
tions, are considered in model development. A representa- 
tive group of the most noteworthy second generation models, 
A ompiled by Grace [16], are listed in Table I. 

The Bubbling Bed Model, developed by Kunii and Levenspiel 
[21] was shown to be semi-successful in correlating experi- | 
mental reaction data and is a useful design tool. However, 
1t utilizes an effective or average bubble diameter to 
describe the bubble phase. This ignores the effect of 
bubble growth on the operation of the bed and, therefore, 
seems το be an oversimplification. “An improvement canepe 
found in the models of Toor and Calderbank [22], Yates et al 
[23], and Wen and co-workers [1,24], where changing bubble 


size is introduced. 
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TABLE I 
SECOND GENERATION FLUIDIZED BED MODELS 


ORG ILI com 
HOVMAND and DAVIDSON [18] 
PARTRIDGE and ROWE [19] 


CHIBA and KOBAYASHI [20] 
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Wen and Kato [24] investigated reactions using their 
Bubble Assemblege Model and found the model predicts conver- 
slons and concentration profiles reasonably well. This same 
model in modified form was used by Wen and Yoshida [1] to 


correlate results of a non-catalytic ore roasting reaction. 
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ΓΙ. MODEL) DEVELOPMENT 


The first step in the development of a model for a 
fluidized bed is the study of the qualitative nature of 
ША таї2ат1їоп. In Figure 2, the three possible statesmor a 
bed of solids having a gas flowing through it are shown. 

At low gas flow rates, Figure 2a, the solids are in a fixed 
Peeecondition. The gas passes through the interstitial 
voids of the solids. As the gas flow is increased, a point 
will be reached where the pressure drop of the gas is suf- 
ficient to support the weight of the bed. At this gas flow 
rate, the bed becomes fluidized and has many of the proper- 
ties Of ordinary liquids. In this state, Figure 2b, the 
bed is characterized by a minimum fluidization velocity, 
voidage, and height. | 

Increasing the gas flow above the minimum fluidization 
rate causes the bed to transition to a state of vigorous and 
vlolent bubbling as shown in Figure 2c. This state is Of 
the greatest importance for it represents the condition of 
most industrial fluid bed reactors. 

In modeling the bubbling fluidized bed, two assumptions 
are universally accepted; 1) isothermal operation and 2) the 
flow through the bubble phase represents all gas flow above 
that required for the minimum fluidization, i.e., (U-UMF). 
These assumptions were used in the model developed i this 
investigation. This model represents a modified form of the 


Bubble Assemblege Model proposed by Wen and Yoshida [1]. 
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A. OVERALL BED MODEL 

The fluidized bed was divided into a number of backmixed 
reactors connected in series as shown in Figure 3. The 
height of each section corresponds to the average bubble 
diameter at that level. The model thus accounts for chang- 
mae bubble diameter within the bed. 

J. Bubble Size and Section Height Calculation 

The bubble diameter was calculated using the rela- 

fon of Kobayashi 125 е ро found that experimental data on 


bubble sizes could be correlated by the equation 


U 
рв = 0.684ррГр MF" | (3) 
ο 


ΙΝ. account for the initial bubble size at the surface of 


the distributor Kobayashi's relation was modified to 
bc 
DB = 0.684ppDp UMF + Do (4) 


where Do is the initial bubble size at the арат 
plate. This quantity is difficult to estimate for beds 
having porous plate distributers and therefore must be 
estimated. For perforated plate distributors, however, Do 
can be estimated based on the work of Davidson and 

Harrison [14]. It was shown by these investigators that 
the bubbles produced in fluidized beds have essentially 

the same properties as bubbles produced in ordinary liquids 
having a small viscosity. On this basis Do can be calcu- 


lated from the equation developed for predicting the 


2 





+ 


Figure 3. Sectionalized Modell or the Fluid (Bed Reactor. 





Solid 
Removal 


EIS 


Solid Feed 


22 





diameter of bubbles produced in ordinary liquids at a 


single orifice 


4 4 
6 — [U-UMF Y 
(=) 


box No g (5) 





The averaging procedure used in calculating the indi- 
vidual section heights is graphically depicted in Figure 
4. For the first section 


Do*XMAH,-*Do 


£ 1 
. 2р 
АН: = (м (7) 


Repeating the procedure for the second section yields 


(2* XM) 


AH 5 
(2-ХМ) 


2 = 2 Do (8) 


In general, it can be shown the average height of the yee 


section is 
АК o (9) 


The bubbles were considered to grow continuously until 
the maximum bubble diameter predicted by Harrison et al 
[26] is reached 


2 
(10) 





= 
Ud 
И 
са | 


This maximum is attained when the upward velocity of the 


bubble equals the terminal velocity of the bed particles. 
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The calculation of the terminal velocity is given in the 
Appendix. 

2. Bed Voidage and Height Calculations 

The bed conditions at minimum fluidization must be 

defined prior to the calculation of the operational bed 
characteristics. These calculations were completed using 
well defined experimental correlations that are given in 
the Appendix. | 

Davidson and Harrison [14] have shown that the rise 
КО без ту of a crowd of bubbles in a fluid bed can be cal- 


culated from the relation 


U,, = U-UMF + 0.711 | (11) 


ie 
DB |^ 
BR 


Š 12 
It was assumed that the expansion of the bed over that at 


Minimum fluidization is a reflection of the total volume 


of bubbles in the bed 
УВт = (Н-НМЕ)5+ (12) 
andthat an average bubble size can be taken as the bubble 


Size at a height equal to HMF/2. The bed expansion ratio 


was calculated by combining equations (4), (11), and (12) 





to give 
VU —— BE (13) 
ШҮ PA 
where 
E HMF 
DB ave = XM F t DO (14) 


u 





This equation was tested by Wen and Yoshida [1] and shown 
to correlate experimental bed heights to within 10%. 
Based on the above result, the bed voidage was calcu- 


lated as 


1-е = a (1-ЄМЕ) (15) 


B. SECTION MODEL 
The sections into which the bed was ο ο were treated 
as having three phases; a bubble phase, a cloud and wake 
phase and an emulsion phase. Figure 5 shows a single bubble 
and the three phases. Flow within each phase was assumed 
to be completely backmixed. This implies uniform gas con- 
centrations and uniform reactant conversion levels within 
each phase. 
ie bubble Phases Description 
The equations used to describe the bubble phase 
are based on experimental results which indicate a bubbling 
fluidized bed can be treated as analogous to a bubbling 
1nviscid liquid. | 
The relative rise velocity of a single bubble was 


calculated from the Davies-Taylor [27] equation 


DB 


5 
U Oz E 23 (16) 


Buc relative rise velocity is the rise velocity of a bubble 
through a bed at minimum fluidization conditions. Ina 
vigorously bubbling bed, the section of the bed ahead of 


the bubble moves upward with a velocity equal to U-UMF. 
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The absolute rise velocity of the bed που Decot: 


U = U = UNE + U 


BR BREL (17) 


This equation was proposed by Nicklin [28] and later by 
Davidson and Harrison [14]. 

An accurate model of any section must predict the 
munber of bubbles within that section. The number of 
bubbles was predicted from the bed expansion ratio according 


to the following logic 


WE Volumemet Section Comprised of Bubbles 
| Volume Per Bubble 


N = Volume ot Section X Bed Expansion Ratio 
Volume Per bubble 


This yields upon substitution 


6 S 
T e - e MF 
N = —— ` —— с (18) 
TAH 2 1-ЕМЕ | 


N 
Maes total volume of bubbles is then given by 


VB = N ' Z рв? (19) 


2. Cloud and Wake Description 
The cloud phase was described according to the model 
of Davidson [29] which permits calculation of the cloud 
diameter according to the relation 


5 EN Όπῃει +20МЕ/ ЕМЕ 


DC” = DB (20) 
Ü.-..- UMF/EMF- 


On this basis, the total volume of cloud phase in each sec- 


tion becomes 
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VC = VB c gU EI ETENE uo 
BREL 

The cloud was assumed to have a voidage equal to the bed at 

minimum fluidization. 

The wake was described using the results of Rowe and 
Partridge [30]. These investigators have shown that the 
wake comprises approximately thirty percent of the bubble 
volume. A plot of their results is given in the Appendix 
and was used. to estimate the parameter Alpha which equals 
the ratio of the wake to bubble volume. The wake was also 
assumed to have a voidage equal to the bed at minimum 
fluidization. 

ο. Emulsmen Description 

The emulsion phase includes all of the remaining 


section volume 


Ms WN S Ye к (22) 


having a voidage equal to the bed at minimum fluidization. 
It was assumed in this investigation that for operating 
conditions where U/UMF is large that the velocity of the 
gas through the emulsion phase is zero. This assumption 
finds substantiation in the work of Latham et al [31] and 
Kunii and Levenspiel [15]. These investigators found that 
the emulsion gas will reverse its direction (Ue-0) when 
Ш ИМЕ > 2.7 = 6.0. 
4. Solid Movement and Gas Exchange 
The primary mechanism of solid circulation within 


the bed is by transport'in the wakes of bubbles. The solid 
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is entrained by rising bubbles and carried upward with a 
velocity equal to the rise velocity of the bubbles. During 
the bubble's movement through the bed, the solids in its 
wake are continuously exchanged with the bulk emulsion 
solids. In this investigation, the exchange coefficient 
for the solids was assumed to be infinite. 

Solids in the bulk emulsion phase move downward 
with the same volumetric flow rate as solids carried upward 
in the wakes. As a result, there is no net flow of solids 
across any horizontal plane within the bed. 

For a system in which a solid is continuousiy fed 
to the bottom of the bed and removed from the top, a net 
upward flow does exist and is equal to the volumetric flow 
rate of the feed. Mathematically, this model can be 


expressed as 


W Qon U 
m BN BRN 
WB = 3, ВМ + 144.0 (23) 
p 
ВМ ВВМ f 
ΜΝ "το a (24) 


where WB(N) represents solids flowing upward to the (КО: 


h 


section from the N** section and WE(N+1) represents solids 


h section from the (N+1 y th 


flowing downward into the nt 
“C Ction. 
The overall gas exchange coefficient was calculated 


using the correlation of Kobayashi et al [32]. Based on a 


unit volume of bubbles, this expression can be expressed as 


Ee] 
Kbe = DB (25) 
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бо :КТМЕТТС5 

Non-catalytic gas-solid reactions in which the reacting 
particles do not change size were studied in this model. 
Examples of this type of reaction include sulfide ore roast- 
[mee reduction of iron ore, land calcination of limestone; 
The following stoichiometric equation was used to represent 


these reactions 
A(gas) + bB(solids) + PRODUCTS (26) 


For this stoichiometric relation, the first order irrever- 


Sible reaction rate equation can be written as 


l dNa _ _ 1 dNb _ aS 
СА A n dt kc Cg (27) 
Yx rx 


The reaction was assumed to proceed according to the 
shrinking core reaction model [1,33]. This model assumes 
the reaction begins as the surface of the material and 
proceeds inward; the reaction zone being the surface of the 
constantly shrinking unreacted core. The particle maintains 
its structural integrity as a result of the formation of a 
completely converted inert ash which surrounds the unreacted 
mone. A time history of the reacting particle is shown in 
Figure 6. It was also assumed that the overall conversion 
rate is controlled by the chemical reaction step. This 
implies that the diffusion of gaseous reactant occurs much 
faster than the reaction at the surface of the core and, 
therefore, a gas concentration gradient between the bulk gas 


phase and the core surface does not exist. 
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p. MATERIAL BALANCE 
Пи. Preliminary Assumptions 
A modified reaction rate constant based on unit 
volume of solids was defined by dividing both sides of 
Equation (27) by Vo to give 
dN 


-1 A u — 
Vo dc Kt CË € 
where 
24k R ° 
Kr = EE (29) 
D 
p 
IE not a constant at a given temperature. It depends 


upon the core radius, Rc, which is a measure of the average 
level of conversion of the solids. It was assumed in this 
investigation that each of the N sections into which the 
bed was divided contains solids having a uniform conversion 
level or stated in terms of the shrinking core model having 
an equal unreacted core size. Furthermore, because of 
backmixing within each section, the solids swept upward 


from the neh 


section to the (N+ 1) th section in the wake 

of bubbles were assumed to have a conversion level charac- 
teristic of the Ln section. This is also true for the 
downward flowing emulsion solids. These assumptions are 
shown in Figure 7. It should be noted that no real physical 
Significance can be attached to the concept of having an 
average core radius within each section. It is obvious 


h 


that the radius of the solids entering the Nt" section from 


the (N+ 1) th section does not increase to that of the material 
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Figure 7. Core Size Distribution. 
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h 


in the N'? section. This assumption was made merely to 


simplify the mathematical representation in terms of the 
section material balances, | 

To further simplify the physical representation, the 
shrinking core model was used to define an average solid 
reactant concentration for any particular solid conversion 
level. This approach can be explained with the aide of 
Fipure 8. It was assumed that within each particle, reac- 
tant is uniformly dispersed through the entire volume of 
the particle. Since the shrinking core model states that 
reaction takes place only at the surface of the unreacted 
core, using the above assumption it can be stated that the 
unreacted core of a partially converted particle has a 
reactant concentration equal to the initial reactant con- 


centration. The average concentration for the entire 


particle was then estimated as 


CS. TC 
CS = — | (30) 
О 
2. Section Material Balances 


Steady state mass balances on the gaseous reactant 
Were used to characterize the flow streams of the system. 
Figure 9 shows an individual section and the associated 
gas flow streams. The material balances for the bubble 


phase and emulsion phase can be written as 


VB [Е CGB] 


EE [CGBy αι – Сов + Ge) N N 


Kr CGB, = 0 (31) 


Y cN N 


5S 
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(о 5) LCCR CE n ЕЕ (32) 


N N 


Included in the reaction term are Ye and να which define 
the distribution of solids between the two phases of the 


section. These terms are defined as 


ха Volume of Solids Dispersed in Clouds and Wakes 
Ша [VC + aVBy] ° [1-=МЕ] /1728.0 (33) 
= Volume of Solids Dispersed in The Emulsion 
ст VEN - [1-eMF1/1728.0 (34) 


The material balances were derived on the basis of zero gas 
flow in the emulsion and complete backmixing within each 
phase. 

in the top reactor section; the cloud phase volume must 
be modified to account for extension of the bubbles above 
the bed surface. The bubble was described as shown in 
Figure 10. in this modification, only the fractional volume 
of the cloud phase associated with that part of the bubble 
below the bed surface is considered. The equations for 


cloud phase then become 


VC = N ° [V Vb 


FR ЕВ! 


where 


DIFH < 0.5 AH, 


= 2 
Vep = 1.0472 [DIFH]“[3.0R.,- DIFH] (55) 
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E | 2 ¿ 
VBap = 140472 1DUEH.R в O DE 
ΙΙ 
_ I ñ 2 а 
УВ р = Мос ОЛГАН pm n) (s Q smo. l h 
E ; я , 2 _ 
VByg VB, -1.0472 [AH- DIFH-R, ,*R5] ^ [2. OR, -AH*DIFH«R, , ] 
(38) 


An overall material balance on both the solid and gas 
Streams was derived from Figure 11 which shows the bed 


th 


representation from the (N-1) section down to the feed 


section. An overall balance yields 


bU [CG -CGB 


vol ди 


ος IWE CS 
О 


N- 11 Е N N 


acc x (59) 


3. Solution of the Material Balance Equations 


A trial and error procedure was used to obtain solu- 
tions to the material balances of the system. The procedure 


can be characterized by the three major steps involved 


Step 1: Assume the exit gas stream concentration CGB(N) 
Step 2: Solve each of the section material balances to 
obtain the required initial gas concentration 
for the assumption of Step 1. 
step 3: Compare the calculated initial gas concentra- 
tion with the known value (CG,) and continue 
the calculation until they are equal by 
adjusting the value of CGB(N). 
Exdetailed description of the completion ot Step 2 199p em 
in the next three sections. 


a. Calculation for the Top Section 


The bed with the yth section isolated is shown 
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Figure 11. Flow Streams for Overall Material Balance. 
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in Figure 12. Four material balances can be derived trom 


this figure: 


U. 1 ÍCGB A 1 -CGBy] + (XK, e) yVBy [CGE 7 CGBy] 
: Y ¿yiryC6By = 0 (40) 
(Ki e) VBy [CGBy "CGE, ] " Y¿yKTyCGEy = 0 (41) 


a Po + MENOS y - NB а ΡΟΝ 


N-1! NCSN Nil 


bU o] [CG ССВ £ 


bU „o1 1EGBy - -CGBy] = WB US A LS. WE CS 


№-1 7 "gESy > WENCSy (45) 


1 Naa 


Equations (40) and (41) are the gas reactant material 

balances for the bubble and emulsion phase, Equation (42) 
is the overall material balance over the πι section 
to the feed section, and Equation (43) is the balance for 


h 


both streams in the nt section. 


The solution of these equations requires knowledge of 


h section for 


the average core size of material in the nt 
evaluationof Kr, the reaction rate constant. This is not 
known and therefore a second trial and error solution pro- 


ceeds as follows. 


Step 1: Guess Rcore and calculate Kr from Equation (29) 
and Vcore. 


Step 2: Calculate CGE(N) and CGB(N-1) from Equations 
(81) апав слоу 


Step 3: —CalculacesCS(N)rfyom Equation 30); 
Step 4: Calculate CS(N-1) from Equation (42). 


Step 5: Repeat the calculations until the section 
material balance Equation (43) is satisfied. 
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Figure 12. Тор Section Material Ба капек 
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When Step 5 has been completed, the solution proceeds to 
the an section. 
b. Calculations for the (Ν-1) 5η Sec lon 
The middle bed sections are characterized by 
the three standard material balance equations which can 


be written as 
Оуо1 СОВу- 27 ССВу.1] + (Къо) м- 178-1 PCGEN_ 4 7 CGBy 1] 


Кт Οδ а=) (44) 


СЕАТ 


С Dis CGE y Pr mco ΠΕΡΙ 


= W CS 


bU o1 [CCo -CCB Ν-15ψΝ-1 


CS E WE 
О 


N-1) £ 


+ н > NO (46) 


Solution of these equations does not require a trial and 
error procedure since the average core size of material in 
the (№-1) th section is fixed by the resultsrof the calcu- 
lations on the cn) th section through 


CS Vcore, 4 


бе 4 = τ (47) 


This equation defines Rcore and Kr for the we)! section. 


solution then proceeds as follows 


Step 5: "GGE(N-T) 5 €GB(N-2) are calculated from 
Equations (45) and (44). 


Step 6: CS(N-2) is calculated from Laudtron (πο. 


This procedure is repeated until the feed section is reached. 
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с.” Calculation “for the Feed Section 
The material balances around the feed section 


can be written as 


U τοι [σος 6681 + (K, 


521 VB, [CGE, -CGB, ] 
- Y¿¡Xr,CGB, = 0 (48) 
O UB CCE ЕЕ hao O (49) 


As is the case for the middle sections, a trial and error 
procedure is not required for their solution on this basis, 


solution was completed as follows 


Step 7: Rcore is calculated from CS(1) obtained from 
the calculations on section two. 


Step 8: Kr and CGE(1) arémealculated from Equations 
(29) and (49). 


Step 9: σα, is calculated from Equation (48). 


The initial gas cos о calculated Step 9 is them 
compared with the known value of σα. This provides a 

test of the validity of the assumed value of the exit gas 
concentration CGB(N). I£ the calculated ος, does not equal 
the experimental value, CGB(N) is adjusted and the entire 
calculation repeated. When the equality exists, the materjial 
balances for each section and the overall material balance 
for the bed have been satisfied and the solution completed. 
The computer flow chart for the solution of the material 


balances is given in Figure 13. 
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Figure 13. Flow Chart for the Solution 
of the Material Balances. 
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КОО THE COMPUTER MODEL 

The program developed to model the fluid bed reactor is 
given in the COMPUTER PROGRAM SECTION. The program was 
designed to be completely self-contained, that is, no 
external calculations are required. The input variables 
are listed in Table II. They represent common design 
variables and therefore would normally be well defined. 
In Table III the program is divided into its three major 


sections and the calculations completed in each are listed. 
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N 


о со N а мл A WN 
e e o% ò o - o 


10, 
12. 
HES 
14. 
πο. 
16. 


Шу, 


13. 
о 
20. 
24. 


ТАВЬЕ 11 


COMPUTER INPUT VARIABLES 


1. Reactor Diameter 


Number of Holes per in.? in the 
Distributor (0 for Porous Plate) 


Weight of Solids in the Bed 
Diameter of the Bed Solids 
Radius of Reactant 

Density of Bed Solids 
Density of Reactant 

Bed Temperature 
stoichiometric Constant 
Kinetic Rate Constant 
Molecular Weight of Solid Reactant 
Solid Feed Rate 

Таста Sedid Concentration 
initial Gas Concentration 


Superficial Gas Velocity at Room 
Temperature 


Superficial Gas Velocity at the Bed 
Temperature | 


Alpha (Volume of Wake/Volume Bubble) 
Initial Section Height 
Initial Guess of CGB(N) 


Increment of Gas Search 
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IV. PRESENTATION AND DISCUSSION OP PD T μι 


The model developed in the previous section was used to 
study a sulfide ore roasting reaction following the stoichio- 


metric equation 


221$ + 30. > 2Ап0 + 250, (50) 


The experimental results obtained by Yagi et al [34] were 
used for evaluating the model. The performance of the model 
was Judged by its ability to predict experimental conver- 
sions. This presentation is divided into two sections. In 
the first, the overall results are presented predicted by 
the model of Wen and Yoshida [1]. In the second section, 
the effects of the major assumptions used in the development 
of the mois are investigated. This section represents a 


detailed study of the model performance. 


'ὶ FOVERALL RESULTS: GENERAL MODEL PERFORMANCE 

In Table IV, the experimental data of Yagi et al [34], 
the calculated conversions obtained by Wen and Yoshida [1], 
and those from the model developed in this investigation 
w tabulated. In these runs, the height of the first 
section AH, was assumed to be 1.0 cm. This assumption was 
proposed by Wen and Yoshida for reactors having porous plate 
distributors. The validity of this assumption is discussed 
in another part of this section. The input data for these 


runs are given in Table VIII of the Appendix. 
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Кип 
Мо. 


10 


11 


12 


15 


Сопу. 


Exp. 


DON 


ΟΠ 2. 


88.7 


91.0 


SOVS 


935 


85.4 


80.6 


72.4 


Sue) 


Сопу. 
Mode1 


НОЕ 


"M 


99.4 


D 


972 


Sho 


одаи 


О? 


81. 6 


оо 


TABLE IV 
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Error 


107 


11. 
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11 
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759 
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25 
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‚01 


62 


‚49 


COMPARISON OF OVERALL RESULTS WITH LITERATURE DATA 


Conv. 


WEN 


ЭОК 


ο 


90.0 


95 0 


δ6. 9 


95.0 


Error 





The data shows that the model predicts conversions which 
are significantly higher than the experimental values. For 
experimental runs having relatively high conversions (>90%), 
the material is essentially completely reacted within the 
bed (>99% conversion). The average percent error of the 
calculated conversions for all runs was 8.3%. The maximum 
error was 12.01% for Run 10 and the minimum was 0.59% in 
Run 1. The model did not predict for any of the runs, a 
conversion less than the experimental value. Low percentage 
error for runs having higher conversion can therefore be 
Ewected. 

The calculated conversions reported by Wen and Yoshida 
show much better agreement with the experimental data. The 
average percentage error for their runs was 1.25$. These 
investigators did not report results for experimental runs 
in which low conversions were found and for SES the model 
proposed in this investigation had the highest error. On 
the basis of the six runs reported by Wen and Yoshida, the 


error for the proposed model becomes 6.60%. 


Ши EFFECTS OF THE MAJOR ASSUMPTIONS 
1. Modified Section Height Calculation 
The height of each section was originally calculated 
by taking the average bubble diameter at a particular height 
as described in the MODEL DEVELOPMENT SECTION. When this 
is done, a portion of the cloud surrounding the bubble 
actually extends beyond the section boundaries as shown in 


Figure 14. The volume of the cloud phase associated with 
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Figure 14. Section Height Based on Bubble Size. 
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each section, therefore, 1s slightly larzer than ТОСО Л sh 
value. Furthermore, this portion of cloud volume is counted 
twice when succeeding sections are considered. 

Figure 15 shows a plot of the gas concentrations for 
the bubble and emulsion phase as a function of height above 
the distributor. It can be seen that the bubble gas concen- 
tration is considerably greater than the concentration of 
gas in the emulsion phase. As a result, solids in the 
bubble wake and cloud react with gas of higher concentration 
than emulsion solids. It was theorized that the inaccura- 
cies in the prediction of the volume of the cloud phase 
resulting from the section height calculation, might con 
(tribute to the prediction of a larger overall conversion. 

To test this theory, the section height calculation was 
modified to include the entire cloud volume. The equations 
necessary to make this modification and their derivation 
Er. given in the Appendix. 

ime conversion results obtained from the modified pro- 
gram are given in Table V. No significant change in the 
predicted conversions was found. The overall атс єл 
was 8.30% which equals that found for the original program. 
The average error for the six runs reported by Wen and 
Yoshida was found to be only slightly higher at 6.64%. The 
maximum change in predicted conversion occurred for Run 12 
where the conversion calculated from the modified program 
was found to be 1.6$ higher than the original calculation. 


These results show that the error in cloud phase volume 
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Run 
No. 


10 


17 


1.2 


13 


А: 
В: 


TABLE V 


PREDICTED CONVERSIONS FOR MODEL WITH 
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mesulting from calculating the section heights on пе a i 
of average bubble size is not a significant factor. In the 
rest of this study, however, the modified section height 
calculation was retained. 

2. Backmixing Considerations 

In the development of this model, it was assumed 

that the solids within a particular section have an equal 
average concentration and conversion level. This concept 
maplicitly assumes that every particle fed to the reactor 
has a residence time which is equal to that of any other 
particle. The mean residence time being defined as 


“OD Fre 51 
a | (51) 


Investigators [15,35] studying the flow patterns and 
residence times in fluidized'beds have shown that a residence 
time distribution (RTD) actually exists for the feed 
particles. This concept implies that some feed solids have 
a residence time which is far less than the average, while 
others have residence times greater than the average. As 
a result, the product of a reaction occurring in a fluidized 
bed does not have a uniform conversion level. Kunii and 
Levenspiel [15] proposed that the RTD can be expressed by 


defining an exit age distribution function E(t) as 


E(t) = emer (52) 
t 


where E(t)dt is the fraction of solids staying in the bed 


tor a time between t and t+*dt. This expression has been 
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shown to correlate experimental residence times very 
accurately. 

The normalized exit age distribution function is plotted 
in Figure 16. The shaded area represents the fraction of 
material having a residence time equal to t+10%. This 
fraction was found to be only 7.2% of the total bed 
material. This observation indicates that the assumption 
that all the solids in the bed have a residence time equal 
to t might be a gross oversimplification leading to poor 
lated results. 

To test this theory, overall bed conversions based on 
a non-uniform bed residence time were calculated using 
average gas concentrations predicted by the model. The 
equations used were those given by Kunil and Levenspiel 
[15). According to these authors, the mean conversion of 
the product stream can be defined as | 

t т | t 
XB = |] E E , 1 


1 


3 са 
1.0 : ot (53) 





where the time for complete conversion is equal to 


E ро 20525) 
ALA | 


The overall average gas concentration Cg was calculated as 


a weighted average of the section gas concentrations. 


N an —— 
L CCE 
a 
een (55) 
t 


The results of the calculation are presented in Table VI. 
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Exit Age Distribution. 
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Кип 
Мо. 


10 


11 


1% 


13 


TABLE VI 


RESULTS OR CONVERSTON CALCULATIONS 
FROM BACKMIXED EQUATIONS 


0, 


οι οσο. 
" Фа Convers. 
99.4 99.8 99.5 0.1 
97.2 98.5 98.9 ΠῚ 
88.7 90.0 95.6 7.78 
91.0 93.0 94.9 4.29 
86.5 86.9 87.8 1.60 
ET 95.0 97.1 3.85 
85.4 š 87.9 — 2.493 
80.6 Е 75.7 6.08 
722 : 61.2 1597 
91.7 - 98.2 7.09 
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The conversions calculated from the backmixed equations 
show considerably better agreement with experimental results. 
The overall average error was reduced to 5.08% and the error 
for the first six runs to 3.21%. These compare with errors 
of 8.30% and 6.64% obtained from the model. These results 
indicate that neglecting the existence of a residence time 
distribution for the bed particles causes the prediction of 
conversions which are high when compared with experimental 
results. 

The Bubble Assemblage Model proposed by Wen and Yoshida 
does account for the non-uniform residence times of the bed 
particles. The close agreement of their results with experi- 
mental data is a further indication that the assumption of 
uniform solid conversion levels in each section incorporated 
in the model proposed mn this investigation is an over- 
simplification. | 

ὄπ Analysis of the Initial™Reactor Section 

One of the major assumptions required in the analysis 
of the model performance was the assumption that the height 
of the first section is equal to 1.0 cm. This assumption is 
necessary because the reactor used by Yagi had a porous plate 
distributor for which an accurate calculation of initial 
bubble size is not possible. The effects of this assumption 
on their model's performance were tested by Wen and Yoshida. 
They concluded that the calculated conversions were not sen- 
се tO Еве Вене of the first section cas tandas ola 
reactant concentration profiles given by these investigators, 


however, show very large changes occurring in the area 
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directly above the distributor. This observation indicates 
that the conclusion that the initial section height does 
NOt effect the overall results might not be adequate for 
all reaction conditions. In this section, the results of 
meectalled study of the first reactor section are given. 

It was found that essentially all the reaction occurs 
1n the first section of the bed. Data supporting this 
observation are given in Table VII. These data show S 
the concentration of solids leaving the first section is 
reduced to very nearly the final product concentration 
level. This observation is not surprising. Many investi- 
gators [36,37] have reported that the primary reaction zone 
is the area directly above the distributor. This fact 
emphasizes the necessity to accurately model the initial 
section. 

Wen and Yoshida have suggested that the only criterion 
required in determining the effects of the first section 
height on the overall model performance is the kinetic 
speed of the reaction, i.e., the value of the rate constant 
kc. They conclude from their results that the overall con- 
version for a fast reaction is not effected by the value of 
АН, and, therefore, no effect will be seen for reactions 
having low rate nsn ns 

In Figures 17 and 18 gas concentration profiles for 
Runs 1 and 12-are plotted. These two runs represent the 
extremes in experimental conversions which were investigated. 
Run 1 is kinetically the faster of the two; having a kc = 


2 2 


ο ο ος ο σος compared With ke = 0.328 x 10 ey ses 
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TABLE VII 
“SOLID CONCENTRATION DATA 


Eun Initial Avg. Conc. Avg. Conc. 
No Conc. 5 Leaving Leaving 
i mole/ft. 1St Section Top Section 

1 NUT 4.3 x 10^ 2.5 χ 105} 
2 στον OE 10 7.6 х 1074 
- 2 - 2 

6 3.58 8.7 x 10 8.6 x 10 
г. 9.9 x 10 2 


10 5. 5898 ια κ 10 
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for Run 12. The figures show that gas concentrations for 
the kinetically slower reaction change more Fons M. in the 
initial height of the bed than do the concentrations for 

the more rapid reaction. The emulsion gas concentration 
showed a very rapid drop for Run 12. At a height of 0.5 in. 
from the distributor, a 71$ decrease in emulsion gas,concen- 
tration was found. This compares to a 27% change in the 
emulsion gas concentration for Run l. 

The effect’ on the overall conversion resulting from 8 
change in the initial section height for Runs 1 and 12 is 
shown in Figure 19. As suggested by the gas concentration 
profiles, Run 12 shows a sensitivity to AH), while Run 1 
does not. This trend is not expected if Wen's and Yoshida's 
conclusion that slow reactions are mot sensitive to the value 
of АН, 1s accepted. 

In Run 12, the reactant feed rate and the initial reac- 
tant concentration of the feed were substantially larger 
than Run 1. The fact that the conversions in Run 12 were 
sensitive to the value of AH, suggests that the feed rate 
and initial concentration as well as the value of the rate 
constant are important parameters in determining the effect 
of AH, on the overall bed conversion. 

Runs 10, 11, and 12 were used to study this effect. 
These runs represent a series in which all parameters 
except the feed rate are constant. The feed rate increases 
for this series of runs. 

In Figure 20, normalized average gas concentrations are 


plotted for the three runs. As the feed rate of the solid 


66 





Overall Conversion 


100. τ 


95. 


90. 


85. 


80. 


Що: 


70. 


65. 


Figure 19. 








0 





0 





Run 1 


Run 12 


. 4 8 1.0 T. ° 


Height of Initial Section (in.) 


Conversion Versus AH; Rugs, 1, 12. 


67 





Normalized Average Gas Concentration 


Eronre 20% 


o Run 10 


^ Run 11 


+ Run ο 





. 4 Ө 
Normalized Height 


Normalized Gas Concentration 
Profiles Κιπσ IT m2. 


68 





to the bed increases, conversion and average gas concentra- 
tions decrease as would be expected. All the runs show 
rapid gas concentration changes near the distributor plate 


which suggests that the value of AH, will effect the overall 


1 
calculated conversion. This was found to be true as shown 
ШГ сиге 21. Кип 12 showed the greatest sensitivity to 

a change in the initial section size. The conversion for 


Run 12 decreased by 21% over the range of AH For Rum 11, 


1: 
the change was only 15% and for Run 10, it was only 5%. 
These results show that the feed rate does effect the sensi- 
tivity of the model to the initial value of АН,. 

The complete meaning of this conclusion is not clear, 

В опе speculation concerning the reasons for this feed 
rate dependency can be made. The fact that conversion was 
шеши to be dependent on the value, o£ AH; can be interpreted 
to mean that the model will predict experimental conversions 
accurately if correct values of АН, are known. This logic 
ignores the fact that the model might still predict incorrect 
conversions due to its not accounting for a residence time 
distribution of solids within the bed, an effect which was 
previously shown to be significant. 

From Figure 21, the values of AH), for which the model 
w111 predict the correct Е АЕ сап Бе 
read. The values are: Run 10, > 1.0; Run 11, 0.88; and 
Run 12, 0.83 in. The interesting fact is that AH, 1s 


different for the three runs and decreases as the feed flow 


rate 1s increased. One would expect that for the same gas 
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flow conditions, bubbles of equal size would form at the 
distributor. The trend, however, suggests that different 
Size bubbles are being formed at the distributor and that 
their size is dependent upon the reaction conditions in 
ine bed. 

Recent investigation [37] into the mechanics of the 
bubble formation at distributors in fluidized beds have 
shown that very small bubbles are initially formed, which 
Peeause of their large inertia swirl in turbulent jets just 
above the distributor surface. These small bubbles rapidly 
coalesce and form large, stable bubbles at the tips of the 
Meets. If this model otf bubble formation is correct, a 
dependence between the size of thebubbles formed at the 
jets and the reaction conditions is not surprising. 

Consider the effects of increasing the feed flow rate 
to the bed. This action increases the tendency for reaction 
and thus the disappearance of gas from the bubble phase. 

As a result, two effects will occur: 1) smaller, stable 
bubbles will be formed by coalesence and, 2) sharper gas 
concentration drops will occur. These trends are both 
present in the data plotted in Figures 20 and 21; a fact 
which supports the speculative logic used in the analysis. 
The conclusion that can be drawn from these results is that 
the uncertainty in the existing mechanism of bubble forma- 
tion could possibly be the controlling factor in the develop, 
ment of a successful mathematical model of fluidized bed 


operations. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The major conclusions of this investigation are: 
1. The conversions predicted by the proposed model are 
significantly higher than the experimental values of Yagi 
et al [34]. The average percent error of the predicted 


conversions was 8.3%. 


2. The results are not dependent on the method used for 
calculating the section heights. Conversions obtained from 
the model when the section heights were based on the cloud 
diameter were essentially identical to the conversions 
predicted by the model when the section heights were based 


on the bubble diameter alone. 


3. The assumption that the ned solids have a uniform resi- 
dence time equal to t, which ignores a solid residence time 
distribution, contributes to the calculation of conversions 
which are higher than experimental values. Calculations 
using the residence time distribution concept and the 
average gas concentration predicted by the model show a 
Significant improvement in the overall results. For these 
calculations, the average percentage error was reduced to 
5.08%. Modification of the model material balances to 


include the residence time distribution concept is recommended. 


4. Investigation into the effects of having to assume the 


, 


height of the initial reactor section indicate that the 
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overall conversion results are dependent upon the assumed 
value (AH;)). Large changes in the solid and gas concen- 
trations in the initial section support this conclusion. 
The sensitivity of the predicted conversions to the value 
of AH, was found to be a function of the kinetic rate con- 
stant and the feed rate of solid reactant. The trends 
indicate that the sensitivity increases for decreasing rate 
constants and for increasing feed flow rates. 

MMe model proposed in this investigation represents 
the first step in the development of a mathematical model 
for use in the design and study of fluidized bed reactors 
for the combustion of solid wastes. Modification to 
include the concept of a solid residence time distribution 
is required. Furthermore, complete understanding of the 
mechanism of bubble formation at the surface of the dis- 
tributor in a fluidized bed is essential in light of the 
observed dependence between the overall conversion and the 
initial section height. For use in the study of combustion 
type reactions, modification for handling the kinetics of 
Shrinking particles must also be made. In this regard, the 
concepts proposed by Kunii and coworkers [15,38] are 


recommended. 
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APPENDIX А 


TERMINAL VELOCITY CALCULATION 


The terminal velocity of the bed particles was calcu- 


lated using the equations derived by Leva [39]. 


2 
(op -0,) 8Dp 


т TETTE Re « 2.0 (56) 
0.152Dpl: 145: 714 (op. AL 
Ша нак 20<Re<500 (57) 
u О 
g 
.5 
δρυρίορ-ρ. ) 
Up = — 4 Re > 500 (58) 
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APPENDIX B 


CHARACTERISTICS AT MINIMUM FLUIDIZATION 


The bed voidage at minimum fluidization was calculated 
from the data of Agarwal and Storrow [41]. The EMF was 
found by these investigators to be a function of the 
particle diameter and the bed material. Their data for soft 
brick particles is shown in Figure 22. This data was fitted 


to three straight lines. The correlation coefficients of 


ии  . The equations are 
Dp < 0.003 
ЄМР = - 50. 00 + 0.615 (59) 


0.003 < Dp < 0.006 


ЕМЕ = 015.8 Dp + 0.542 (60) 
Dp > 0.006 
MF = -4.18 Dp + 0.474 | (61) 


The minimum fluidization velocity was calculated from 
an equation given by Kunii and Levenspiel [15]. For small 


particles this equation is 


2 
Dp (0Pp-0,)8 | 
UMF = —Ə—ns (62) 
1650μ' 


The gas viscosity and density required for this calculation 
were calculated at the bed temperature and atmospheric 


pressure from the equations 


og = ua IDEAL GAS LAW (63) 
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is 3/2 τς. 
u шш d (64) 


D SE 


1; 
y (М. ) Q 


s (65) 
m1x » y; (M. ) ° 


Equations (64) and (65) were obtained from Perry's Handbook 
for Chemical Engineers [41]. 
The height of the bed at minimum fluidization was cal- 


culated from the EMF by the equation 


12 Ho 
Ec 1.0-eMF > 
where 
144 W. 
Ho = Ost (67) 


f 





Vw/ (Vw+V,) = a/(1+0) 


e 
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Figure 23. Plot of Data for Determination 
of the Wake Volume. 
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APPENDIX D 


TABLE VIII 
INPUT DATA 


We Co 
x 10° μ΄ ΠΙΟ; 
1b/sec mole/ft 


CS, 


3 


πο... 


0.05 
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APPENDIX E 
MODIFIED SECTION HEIGHT CALCULATION 


The cloud diameter surrounding a bubble of size DB can 


pesca lculated from the relation proposed by Davidson [29] 


Upp + 2 UMF/eMF]'/ > 
ВОС ие (68) 
BR 
In order to define a section height (AH) on the basis of 
the cloud diameter, an averaging procedure similar to that 


employed previously was used. Using this procedure, the 


average section bubble diameter was defined as 


N 
X > H + Do 
DB, - oA. ON ` (69) 
Ет 


On this basis, the section height based on the cloud 


diameter can be written as 


5 5 € 
(70) 


U 
_ BR 
AHy = DBy RIAM 


When this formulation is used, the number of bubbles in 


each section becomes 


6 S.. (H-HMF) AH 
N = N (71) 


TH DBy 
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